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COMPUTATION OF BRANCH METRIC VALUES IN A DATA DETECTOR 

FIELD OF THE INVENTION 
The present invention relates to data detectors. More particularly, the 

5 present invention relates to the computation of branch metric values in a data 
detector, such as Viterbi-like detector employed in a read channel of a disc drive. 

BACKGROUND OF THE INVENTION 
A typical disc drive includes one or more magnetic discs mounted for 
rotation on a hub or spindle. A typical disc drive also includes a transducer 

10 supported by a hydrodynamic air bearing which flies above each magnetic disc. 
The transducer and the hydrodynamic air bearing are collectively referred to as a 
data head. A drive controller is conventionally used for controlling the disc drive 
based on commands received from a host system. The drive controller controls 
the disc drive to retrieve information from the magnetic discs and to store 

15 information on the magnetic discs. 

An electromechanical actuator operates within a negative feedback, closed- 
loop servo system. The actuator moves the data head radially over the disc 
surface for track seek operations and holds the transducer directly over a track on 
the disc surface for track following operations. 

20 Information is typically stored in concentric tracks on the surface of 

magnetic discs by providing a write signal to the data head to encode flux 
reversals on the surface of the magnetic disc representing the data to be stored. In 
retrieving data from the disc, the drive controller controls the electromechanical 
actuator so that the data head flies above the magnetic disc, sensing the flux 

25 reversals on the magnetic disc, and generating a read signal based on those flux 
reversals. The read signal is typically conditioned and then decoded by the drive 
controller to recover data represented by flux reversals stored on the magnetic 
disc. 

A typical disc drive read channel includes the data head, preconditioning 
30 logic (such as preamplification circuitry and filtering circuitry), a data detection 



STL 11305 

-2- 

and recovery circuit, and error detection and correction circuitry. The read 
channel can be implemented either as discrete circuitry, or in a drive controller 
associated with the disc drive. 

Some conventional disc drive read channels employ data detection 
5 schemes that are designed under the assumption that additive white Gaussian 
noise is present in disc drives. However, it has been observed that media noise in 
disc drives is neither white nor stationary. The non-stationarity of the media noise 
results from its signal-dependent (or data-dependent) nature. Consequently, 
some more recently developed data detection schemes employed in read channels 
10 have utilized data-dependent noise prediction to account for the inherent data- 
dependence of media noise. However, these schemes suffer from the burden of 
requiring a large number of parameters to be estimated or tuned within the read 
channel. 

Embodiments of the present invention provide solutions to these and other 
15 problems, and offer other advantages over the prior art. 

SUMMARY OF THE INVENTION 
The present embodiments relate to a data detector in which branch metric 
values are computed using a transition jitter model (dependent upon positions of 
data transitions) of media noise, which results in a reduction in a number, of 
20 parameters to be estimated in the detector, thereby addressing the above- 
mentioned problems. 

A method of determining branch metric values in a detector is provided. 
The method includes receiving time variant signal samples, and computing the 
branch metric values as a function of transition jitter statistics corresponding to 
25 the signal samples. A detector configured to determine branch metric values as a 
function of transition jitter statistics corresponding to signal samples is also 
provided. 
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Other features and benefits that characterize embodiments of the present 
invention will be apparent upon reading the following detailed description and 
review of the associated drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is an isometric view of a disc drive. 

FIG. 2-1 is a simplified block diagram of a read channel of the disc drive 
shown in FIG. 1. 

FIG. 2-2 is a block diagram of a data detection and recovery circuit 
according to the present invention. 

FIG. 3 is a simplified block diagram illustrating a generic branch metric 
calculation unit in accordance with an embodiment of the present invention. 

FIG. 4 is a conceptual flowchart showing branch metric computation for a 
general state transition in accordance with an embodiment of the present 
invention. 

FIG. 5 is an exemplary trellis diagram for illustrating an embodiment of the 
present invention. 

FIG. 6 is a table including decimal state representation information for 
illustrating an embodiment of the present invention. 

FIGS. 7 and 8 are plots illustrating a comparison of results obtained using 
prior art branch metric computation techniques and branch metric computation 
techniques of the present invention. 

DETAILED DESCRIPTION OF ILLUSTRATIVE EMBODIMENTS 

FIG. 1 is an isometric view of a disc drive 100 in which embodiments of the 
present invention are useful. The same reference numerals are used in the various 
figures to represent the same or similar elements. Disc drive 100 includes a housing 
with a base 102 and a top cover (not shown). Disc drive 100 further includes a disc 
pack 106, which is mounted on a spindle motor (not shown) by a disc clamp 108. 
Disc pack 106 includes a plurality of individual discs, which are mounted for co- 
rotation about central axis 109. Each disc surface has an associated disc head slider 
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110 which is mounted to disc drive 100 for communication with the disc surface. 
Surfaces of disc 106 are usually divided into zones, with each zone including 
multiple adjacent tracks. In the example shown in FIG. 1, sliders 110 are supported 
by suspensions 112 which are in turn attached to track accessing arms 114 of an 
5 actuator 116. The actuator shown in FIG. 1 is of the type known as a rotary moving 
coil actuator and includes a voice coil motor (VCM), shown generally at 118. Voice 
coil motor 118 rotates actuator 116 with its attached heads 110 about a pivot shaft 
120 to position heads 110 over a desired data track along an arcuate path 122 
between a disc inner diameter 124 and a disc outer diameter 126. Voice coil motor 
10 118 is driven by servo electronics, which is included in control circuitry (or 

controller) 130, based on signals generated by heads 110 and a host computer (not 
shown). 

FIG. 2-1 is a simplified block diagram of a read channel 200 of disc drive 
100. For simplification, only one disc and one head are shown in FIG. 2. Read 
15 channel 200 includes magnetic disc 106, data head 110, preconditioning logic 202, 
data detection and recovery circuit 204 and error detection and correction circuit 
206. Preconditioning logic 202, data detection and recovery circuit 204 and error 
detector and correction circuit 206 are, in some embodiments, a part of control 
circuitry 130 (FIG.l). 

20 As mentioned above, in operation, controller 130 receives a command 

signal from the host system which indicates that a certain portion of disc 106 is to 
be accessed. In response to the command signal, servo electronics within 
controller 130 produces control signals that direct voice coil motor 118 to rotate 
actuator 116 and thereby position head 110 over a desired track. 

25 Head 110 develops a read signal indicative of flux reversals in the track 

over which head 110 is positioned. The read signal is provided to preconditioning 
logic 202 which typically includes a preamplifier, an analog to digital converter 
and filtering circuitry. The amplified and filtered signal is provided to data 
detection and recovery circuitry 204 which recovers data encoded on the surface 
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of disc 106. Once the data is detected and recovered, it is provided to error 
detection and correction circuitry 206 which may be based on an error correction 
code (ECC), such as a Reed-Solomon code. Error detection and correction circuit 
206 detects whether any errors have occurred in the data read back from the disc. 
5 Further, in some embodiments, error detection and correction circuit 206 is 

provided with error correction logic which is used to correct errors discovered in 
the data read back from disc 106. The corrected data is provided to the host 
system. 

Data detection and recovery circuitry 204 typically includes a data detector 

10 (such as a Viterbi-like detector) that helps recover data from the readback signal. 
Operation of a Viterbi-like detector is more easily understood using a trellis 
diagram (such as the trellis diagram shown in FIG. 5), which is a typical state 
machine diagram with an additional parameter, discrete time. The Viterbi-like 
detector operates by selecting a most likely path through the trellis diagram given 

15 some received sequence. A "metric" is kept for each state at each time, and a 
"previous state" is also kept for each state at each time. As new samples are 
received, new metrics are computed. 

As mentioned above, some conventional disc drive read channels employ 
data detection schemes that are designed under the assumption that additive 

20 white Gaussian noise is present in disc drives. In such schemes, trellis/ tree 

branches are usually computed as Euclidian metrics. In general, the bit error rate 
(BER) performance of a detector employing a Euclidian metric computation 
method is relatively low. As noted above, certain other prior art data detection 
schemes utilize data-dependent noise prediction to account for the inherent data- 

25 dependence of media noise. However, these data detection schemes suffer from 
the burden of requiring a large number of parameters to be estimated or tuned 
within the read channel. 

Under the present invention, a scheme for determining branch metric 
values in a detector is provided in which branch metric values are computed 
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using a transition jitter model of media noise. This results in a reduction in a 
number of parameters to be estimated in the detector. The present invention 
differs from previous solutions in that noise statistics (related to amplitude- 
distortion of signals received by the detector) are not explicitly estimated for each 

5 hypothesized data sequence corresponding to a particular trellis branch. Instead, 
a maximum of only two parameters need be estimated for a particular head and 
zone combination. These two parameters, along with a hypothesized data 
sequence, and an equalized transition response uniquely determine all required 
branch metrics, and implicitly determine the overall noise statistics (transition 

10 jitter noise and amplitude-related noise) corresponding to each branch. In 

addition, the metric naturally exploits the non-causal characteristics of transition 
jitter. The calculation of parameters and the determination of branch metric values 
in accordance with the present invention are described further below. The method 
of determining branch metric values in accordance with the present invention can 

15 be used to provide soft or hard decisions, in both trellis and post-processor 

architectures. A data detector, which implements the branch metric computation 
scheme of the present invention, is described below in connection with FIG. 2-2. 

FIG. 2-2 shows a block diagram of data detection and recovery circuit 204 
in accordance with an embodiment of the present invention. While data detection 

20 and recovery circuit 204 will typically include conventional pulse detection and 
qualification circuitry, it also includes finite impulse response (FIR) filter 208 and 
Viterbi-like detector 210. In some embodiments, in addition to a primary detector 
(such as Viterbi-like 210), a post processor 212 is included to refine the output of 
the primary detector. In designing Viterbi-like detector 210 of the present 

25 invention, effects of wide-band additive noise and media jitter in samples input 
into detector 210 are taken into consideration. An example algorithm suitable for 
implementation in Viterbi-like detector 210 and/ or post processor 212 is 
described below in connection with Equations 1-38. 
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In the example algorithm, the derivation of jitter-noise metrics is carried 
out from a Bayesian viewpoint, where transition jitter is treated as a random, 
nonlinear, nuisance parameter. The example algorithm is described below by first 
developing an appropriate background and model notation. This is followed by a 

5 general discussion of the proposed Bayesian approach and the simplified first- 
order Taylor series model for jitter. An optimal Bayes cost function for first-order 
jitter, which is inherently non-recursive (not implementable in a trellis search 
structure), is then derived. This is followed by the derivation of a recursive branch 
metric that can be modified, as discussed in connection with FIGS. 3-6, for 

10 practical implementation in a Viterbi-like detector. 

1) Discrete Time Media Jitter Model 

Let ak e {+/- 1 } represent a non-return to zero (NRZ), bipolar sequence of 
encoded, precoded data to be written to the disk, with the transition sequence bk 
15 defined by 

Equation 1 

Assume a symbol-rate sampled, discrete-time equivalent, equalized, transition 
response is given by the sequence gk, so that the combined effects of wide-band 
20 additive noise and media jitter can be modeled at the input to a detector as the 
received sample 

Equation 2 

where yi is the normalized, random jitter parameter associated with the I th 
25 transition symbol bi, and rik represents the contribution of wide-band, additive, 
Gaussian noise with variance a] . Without loss of generality it can be assumed 
that rik is white, by simply assuming that n has been effectively filtered so as to 
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remove any correlation in tik, and thus all equalization is assumed to be absorbed 
in gjt. Note that for practical direct current (DC)-free inter-symbol interference (ISI) 
channels (e.g. longitudinal recording) gk will have a finite support. However, it is 
shown later in the application how the proposed Bayesian cost function is easily 
5 modified for the case of an arbitrary ISI channel that does not have a DC null. 

2) Bayesian Marginal Approach for First Order Jitter 

Let the column vector r = [ro, n, . . . , yma] 1 represent M received samples of n, due 
to the transmission of the N <M transition symbols b = [bo, b\,. . . , bti-i ] T . 
10 Corresponding to the transition symbol vector, is the (considered random) vector 
of jitter parameters y = [yo, yi, . . . , yN-i] T - The likelihood function, conditioned also 
on y, is thus 

p{r\b, T, a\ ) = n p(r k \b, y, a 2 „ ) 

Equation 3 

15 If the Bayesian viewpoint of treating ^as a random nuisance parameter with the 
assumed prior density function p(y/b) is adopted, the choice is either to consider 
the joint maximum a-posterori (MAP) estimates of b and y, or to integrate out the 
nuisance parameter from Equation 3, resulting in the marginal: 

p(r\b,(T 2 n )= \ y p{r\b,y,cjl)p{y\b)dy 

20 Equation 4 

Here, the marginal approach is chosen because it requires estimation of the 
transition symbol sequence only, whereas the joint MAP estimation of b and y 
essentially requires the estimation of two parameters, one of which is non-linear, 
for each observed sample. 

25 Given the preference for a marginal approach, however, leads to the need 

for dealing with a major obstacle. Note that g is a non-linear function of y, and 
therefore a closed form expression for Equation 4 does not exist in general. This 
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situation is remedied by utilizing a first-order Taylor series jitter model to replace 
Equation 2: 



n-\ 



/=0 



Sk-i Yi Sk-i 



Equation 5 

5 where 

Equation 6 

where def represents 'defines/ For notational convenience, the linear part of 
Equation 5 can be written as 



10 ^ b lgk-l=gl b 
arid the jitter contribution is expressed as 

zbirig k -!=r T Ck(b) 



Equation 7 



Equation 8 



15 where 



Sk def [gk>g k -\>'">gk-N+\\ 



Equation 9A 



g k def 



g k> § g k-N+\ 



Equation 9B 



20 and 



c k {b)d 5 fb&g k 



Equation 10 



where ' 0 represents vector element-by-element multiplication. 
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3) Non-Recursive Solution 

Assuming that each transition is associated with an independent, and identically 
Gaussian-distributed, y yields 

p(r\ b ) = pir) 

N 

= \1k) 2 cr N exp 



V J 

Equation 11 

With Equation 11 and the uncorrected, Gaussian-distributed, assumption for n/t, 
Equation 4 can be solved and the following cost function can be derived 

v*-Q>) dg^ -log (p(r\b, a*]) 

= a 2 n \og\R(b)\ 

+ (r-Gb) T \R(b)\~\r-Gb) 



Equation 12 



where 



4b) def I K 



C(bfc(b) 



C{b) def [coC,...^.,], 



Equation 13 
Equation 14 



and 
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G def 



g 



S M -1 

Equation 15 

Note that the cost function given by Equation 12 is not recursive because the 
inverse of the data-dependent correlation matrix R(b) is in general not diagonal, 
5 except for the trivial cases of b = 0 or <j\ = 0. Also, for a] = 0, Equation 12 reduces 
to a conventional or standard Euclidean metric 

V(b)=r-Gb\ 2 

Equation 16 

One method for converting Equation 12 into a recursive form is to approximate 
10 R(b)- J with an L-banded matrix, meaning that it has 2L + 1 non-zero 

diagonals. This is equivalent to assuming that media and electronic noise 
contributions are lumped together as a data-dependent, L "'-order autoregressive 
(AR) process, and results in corresponding data-dependent noise prediction 
architecture with L noise-predictive taps. The downside to this approach is that 
15 there is no straightforward method for connecting the large number of data- 
dependent parameters required for detection to known quantities a-priori, and 
thus requires the estimation of these parameters for every head and zone 
combination. 

However, it is noted that R(b) is determined completely by the data 
20 transition sequence b (the multi-dimensional vector with which the cost function 
is searched over), the a-priori known transition response derivative gk (via 
Equations 14 and 10), errand <j] . The jitter variance parameter a] can be deter- 
mined relatively accurately a-priori for a given magnetic media formulation and 
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zone BPI (bits per inch), and thus a Bayesian detector may require only the 
tuning/ estimation of a\ . 

This observation serves as motivation to modify the Bayesian approach to 
derive a recursive media-noise cost function that requires only the knowledge of 
5 o^andcr 2 . 

4) Recursive Solution 

One approach to deriving a recursive cost-function is to neglect all off -diagonal 
terms in Equation 13. This is equivalent to using the following as approximations 
1 0 for Equations 3 and 4 

:P(r k \ b,a 2 n ) = Jp(^|6 3 r^ w 2 )p(r|^Vr 

r 

Equation 17 

p{r\b,a 2 n ) «Ylp(r k \b,ol) 

k 

Equation 18 

15 Solving Equation 17 results directly in the branch metric (to be minimized) 

def 



Mr k9 b)- -\ogi(r k \b 9 o*) 



Equation 19 

Equation 19 is further simplified by ignoring the leading bias term, thus resulting 
20 in 

n 

Equation 20 
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An attractive feature of Equation 20 is that it requires knowledge only of the ratio 
of jitter to additive noise variances, and therefore this cost function may be 
interpreted as a branch metric that can be made increasingly robust to transition 
jitter by simply increasing this single parameter. 

5) Implementation as a Modified Viterbi Detector 

Equation 19 is examined in more detail for practical implementation in a Viterbi 
detector. In particular, by accounting for the finite support of gk and gk, the 
dependence of A(.) on the entire block of N symbols b can be dropped. It is first 
assumed that the equalized transition response g is causal, with I nonzero terms 
(later, this assumption is relaxed and the unit pulse response for DC-content 
channels is used). The derivative of the transition response will, in general, 
contain both causal and non-causal terms. A new, length (h + h + l) transition 

derivative response vector is defined: 

• • • • • • 

gdef[g_ Ix ^g_ v g^g v .^g [2 ] T 

Equation 21 

and assuming that h>l-l new length vectors (h + h + 1) for the transition 
sequence and transition response are defined as 

b k dg^ [b k+Jx ,.~, b k+}9 b k _ l9 ..., b k _ h Y 

Equation 22 

and 

g d<£ [0 V ..,0, g 0 , g, g I _ ] ,0 V -,0f 

Equation 23 

where the transition response vector has h leading zeros before the go term, and h 
- 1 + 2 zeros following the gi-i term. 

With the vectors defined above, the branch metrics given by Equations 19 
and 20 are re-written in the, more appropriate, recursive forms as 
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Equation 24 



and 



A (r k ,b k ) = — — 



5 Equation 25 

where 

c k d£b k Qg 

Equation 26 

As mentioned above, the equalized transition response vector g may not 
10 actually have finite support, in which case, with the (assumed finite support) pulse 
response pi and the NRZ data sequence ak, the noiseless the ideal channel output 
term in Equation 24 is computed as 



big = Y^sA-i 

i-o 
1 7-1 



2 /=o 
/ 

/-o 

Equation 27 

Likewise, the remaining terms of Equation 24 are converted to the NRZ domain 
by writing 

1 h .2 

4 /=-/, 

20 Equation 28 



STL 11305 

-15- 

Note that the inner product in Equation 28 represents the norm-square of cu, and 
thus satisfies 

c[c k >0 

Equation 29 

5 Also, note that this term is exactly zero only for the case of no transitions: 

2 

bk = 0, and thus the overall inverse weighting term (1 + (- J j)c k c k ) imposes a 

higher penalty in Equations 24 and 25 for the case of either an increasing <j 2 Y , 
and/ or for more transitions in bk. Conversely, this weighting term goes to unity if 
either a 1 = 0, or bk = 0, in which case, both Equation 24, and Equation 25 become 
10 the standard Euclidean metric for an ISI channel with no transition jitter: 

A(r k ,b k )Euc = (r k -b T h gf 

Equation 30 

An examination of Equation 28 suggests the following definition of a trellis 
detector state at time k: 

15 Zk ~ ( a k-I l -\? a k-1 2 i-''? a ki a K+\i a ic+Ii-\> a k+I x > ) 

Equation 31 

and shows that Equations 24 and 25 require a trellis with 2<h + I 2 + 1 > states, an ex- 
pansion over the 2 1 states required if a 2 y = 0. Note also that the non-causal effects 

of transition jitter are reflected in Equation 31. Corresponding to the above state 
20 definition, there are 2<h +I 2 + 2) possible transitions between states at time k: 

T k = (x k9 x k+l ) 

~ ( a k-I 2 -\ ' a k-I 2 Q k 9 a k+\ ' a k+l y -\ 9 a k+I l ) 

Equation 32 

and thus the branch metrics in Equations 24 and 25 can be explicitly denoted as 
25 functions of the current equalized sample r*, and the state transition Tk as, 
respectively, A(rk,Tk) and A (n,Tk). FIG. 3 illustrates the generic branch metric 



STL 11305 



-16- 



calculation unit, and FIG. 4 shows a conceptual flowchart that further describes 
the branch metric computation for a general state transition. Inclusion of dashed 
path 402 in FIG. 4 results in the branch metric given in Equation 24. If dashed path 
402 is ignored, the result is Equation 25. 

5 

6) 8-State Modified Viterbi Example 

As in a conventional binary-symbol Viterbi trellis search algorithm, there 
are two branch metrics going into each trellis state, reflecting the local likelihood 
of one unique state transition over another, and these branch metrics are summed 
10 over the length of the trellis to represent an overall path metric. For this example, 
an 8-state trellis (shown in FIG. 5), with 2 = 2, Jj = 2, and h = 2. The decimal 
number associated with each state in FIG. 5 is related to the state definition in 
Equation 31 by the table shown in FIG. 6. 

From Equation 27 the ideal, linear, equalized channel output can be 
15 represented as a function of state transition Tk=(ajt-2, ak-i, cik+i) as 

? 1 

b k g = - [g 0 K - « JM ) + #1 K-l ^ <*k-2 )1 



Equation 33 

and from Equation 28> the norm-squared term c[c k as 

1 * 2 * 2 * 2 

c/'c* =-[^,(%i ~ a kf +SoK -**-i) 2 ~^- 2 ) 2 ] 

20 Equation 34 

At the top of FIG. 5 the two possible transitions from time k that can lead to state 
0 at time k + /, denoted here for simplicity as Ak (0, 0), and Ak (4, 0), are shown. 
Note that in this setting, Ak(ij) signifies the branch metric resulting from a 
transition at time k from state i, to state; at time k + 1. From the table in FIG. 6, Ak 

25 (0, 0) can be identified with the transition in NRZ Tk = (-1, -1, -2, -2), and Ak (4, 4) 
with can be identified with the transition Tk = (+2, -2, -2, -1). Thus Equations 33 or 
34 can be used in Equations 24 or 25 to obtain 
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A*(0,0) = rl 



K k (A$) = g.Mogfl + ^g.V (r * + f l)2 



Equation 35 



a+K)g.) 



Equation 36 



or 

A* (0,0) = rl 

(r k +g t ) 



Equation 37 



2 



A, (4,0) = ■ 

a+(-r)gi 

Equation 38 

10 7) Performance Comparison 

FIGS. 7 and 8 are plots of results obtained from longitudinal recording 
simulations with a Lorentzian transition response model. Gaussian-distributed 
jitter parameters are generated for each transition to represent media noise. In 
obtaining the overall channel signal-to-noise ratio (SNR), the total received noise 

15 power used is the sum of discrete-time signal variances due to additive white 
Gaussian noise, and Gaussian transition jitter, as observed through an ideal low- 
pass filter (transition width variation is ignored here). In FIGS. 6 and 7, the 
vertical axis represents BER and the horizontal axis represents SNR in decibels 
(dB). Plots of FIGS. 6 and 7 represented by solid lines correspond to the 

20 conventional Euclidian metric and plots represented by dashed lines correspond 
to the recursive Bayes technique of the present invention. Specifically, FIG. 7 
compares the BER performance of the Euclidean metric (Equation 30) versus the 
recursive Bayesian metric (Equation 25) for a symbol density of 2.25 and a 
jitter/ electronic noise mix of 50/50. A gain of about 0.5 dB is demonstrated for the 

25 recursive Bayesian metric over the Euclidean metric at about 1E-5 error rates. FIG. 
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8 shows BER results where the mix of jitter/ electronic noise is increased to 80/ 20. 
Here the recursive Bayesian metric demonstrates a gain of about 1.25 dB at le-5. 
The results obtained from the longitudinal recording simulations show that the 
BER performance of a detector employing the recursive Bayesian metric 

5 computation technique of the present invention is substantially better than the 
BER performance of a detector employing the prior art Euclidian metric 
computation method. 

It is to be understood that even though numerous characteristics and 
advantages of various embodiments of the invention have been set forth in the 

10 foregoing description, together with details of the structure and function of 
various embodiments of the invention, this disclosure is illustrative only, and 
changes may be made in detail, especially in matters of structure and 
arrangement of parts within the principles of the present invention to the full 
extent indicated by the broad general meaning of the terms in which the 

15 appended claims are expressed. For example, the particular elements may vary 
depending on the particular application for the data detector while maintaining 
substantially the same functionality without departing from the scope and spirit 
of the present invention. In addition, although the preferred embodiment 
described herein is directed to a data detector for a read channel of a disc drive 

20 data storage system, it will be appreciated by those skilled in the art that the 

teachings of the present invention can be applied to data detectors employed in 
other systems, without departing from the scope and spirit of the present 
invention. Further, the data detector may be implemented in hardware or 
software. The disc drive can be based upon magnetic, optical, or other storage 

25 technologies and may or may not employ a flying slider. As mentioned above, 
transition jitter is a relatively dominant component of media noise and is 
dependent upon positions of data transitions. Transition jitter statistics include 
statistical data corresponding to transition jitter, such as transition jitter variance, 
used in the above equations. 



